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We use lattice simulations to examine the detailed dynamics of inflaton fragmentation during 
and after preheating in X(f> 4 chaotic inflation. The dynamics are qualitatively similar to preheating 
after m 2 </> 2 inflation, involving the exponential growth and subsequent expansion and collision of 
bubble-like inhomogeneities of the inflaton and other scalar fields. During this stage fluctuations of 
the fields become strongly non-Gaussian. In the quartic theory, the conformal nature of the theory 
allows us to extend our simulations to much greater times than is possible for the quadratic model. 
With these longer simulations we have been able to determine the time scale on which Gaussianity 
is restored, which occurs after a time on the order of a thousand inflaton oscillations. 
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After inflation the energy density of the nearly homogeneous inflaton field decays into fluctuations of the inflaton 
and other fields. In many simple models of inflation this process begins with an exponentially rapid stage of decay 
\Q ' that produces highly inhomogeneous, nonthermal field fluctuations. In large field chaotic inflation models preheating 
7-H , occurs via parametric resonance, which has been analyzed both analytically (see e.g. and numerically H • This 
preheating stage is followed by a short, violent transition that leads to a regime of Kolmogorov turbulence Q. 

In 4] this transition stage was studied for a model with inflaton potential V = (l/2)m 2 2 . It was found there 
that parametric resonance leads to growth of fluctuations in the peaks of the initial random gaussian field, giving 
rise to a quasi-stable standing wave pattern of bubbles and nodes. This growth continues until backreaction makes 
I . parametric resonance inefficient, after which these bubbles expand and collide, thus bringing the entire space into the 
y— ( ' strongly inhomogeneous regime. It was found there that the fluctuations produced during preheating are strongly 
non-Gaussian, and that this non-Gaussianity persists long after the end of parametric resonance. 

For technical reasons discussed below, the simulations performed in [4J could not be continued long enough to 
determine the ultimate fate of these non-Gaussian perturbations. In this paper we explore these same questions in a 
■ quartic model with an inflaton potential V = (l/4)A</> 4 . In this model it is possible to run the simulations to much 
ij later times and we were thus able to see the field statistics return to Gaussianity. 

In section|TT]we describe the model we are using and our numerical simulations. In section Hill we present the results 
' of our simulations and discuss their implications. 
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II. THE MODEL 

We consider the potential 

V = \\^ + \gW (1) 

where 4> is the inflaton and \ is another scalar field that is coupled to it. We use LATTICEEASY [| to solve the 
classical equations of motion for the two fields 

4> + 3Hc\ - lv 2 + (A0 2 + .gV) 0-0 (2) 

X + 3Hx-^V 2 x + g 2 (b 2 X = 0. (3) 

Oscillations of the zeromode of <f> after inflation lead to parametrically resonant amplification of modes of x within 
certain resonance bands. These amplified modes of x m turn excite fluctuations of <f>. 

LATTICEEASY uses a comoving grid, meaning the wavelength of produced fluctuations remains constant in pro- 
gram coordinates. In a quadratic model this comoving growth causes a problem for compuater simulations because 
the physical wavelength m _1 at which modes are preferentially produced shrinks in these comoving variables. As 
the universe expands the grid spacing needs to remain small enough to accurately capture this physical length scale. 
Once the universe has expanded enough to violate this criterion, it is no longer possible to continue the simulation 
with any accuracy. 
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In the quartic model, however, the parameters A and g are unitless and there is no fixed physical length scale in 
the problem. By defining a new set of variables 4> pr = 4>/a and r = at the equations of motion can be recast as 



(t>" pr - V 2 (j) pr + (A0 2 r + g 2 x 2 pr ) 4> pr + A = 

A = 



x" 



V Xpr 



2/2 
9 (PprX-pr 



(4) 
(5) 



where primes represent differentiation with respect to r and A represents derivatives of the scale factor that vanish 
for a radiation equation of state, an approximation that is very accurate for this model. [9j In other words we can 
scale expansion out of the equations, and can thus simulate to much later times than would be possible for a quadratic 
model. 

The simulations shown in this paper were done on a three dimensional grid of 128 3 or 256 3 gridpoints. The 
simulations start at the end of inflation when the mean value of the inflaton is </>o = 0.342Af p . Times are reported 
in units of l/(\/A0o) and field values are reported in units of 4>q. All of the results shown in this paper are for 
A = 9 x 1CT 14 and g 2 /\ = 200. 

III. INFLATON FRAGMENTATION: RESULTS AND CONCLUSIONS 

Figures [T] and [2] show the growth of fluctuations of the fields cf> and %■ The information in these plots is well known. 
The field x grows exponentially. Some time after the growth of x starts fluctuations of the field <j) begins growing 
with twice the exponent of x- These plots are included here primarily as a reference for seeing where in the process 
of preheating the fields are at each of the times shown in the plots below. 
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FIG. 1: Variances of the fields <j> (red) and \ (blue). 
(For black and white viewing, \ begins growing ear- 
lier; 4> begins later but quickly catches up to x ) 



FIG. 2: Number density for <f> (red, solid) and x (blue, 
dashed) . 



Figure [3] shows the growth of fluctuations of the x held during preheating. For clarity the figures show values 
only on a one dimensional slice through the lattice. The initial conditions used for the model are gaussian random 
fluctuations with expectation values 



(\4>k\*) = 



l 



2jk 



(6) 



l eff 



intended to simulate quantum vacuum fluctuations See the LATTICEEASY documentation for more details. 

As parametric resonance begins fluctuations of the field begin to grow exponentially. Note that the vertical scale on 
the different frames of [3] is not constant. Since parametric resonance only excites modes with momenta below a certain 
cutoff (fc* ~ (<7 2 A) 1 / 4 ^)o the short wavelength fluctuations rapidly become insignificant. The longer wavelength 
modes, however, remain almost entirely unchanged except for their overall scale. In other words the spatial distribution 
of the fluctuations produced during preheating simply mimics the spatial distribution of the infrared modes of x that 
were present before preheating. 

Fluctuations of <fi are not produced until later, but when they appear they grow due to interactions with the amplified 
X fluctuations, so the peaks of 5<p mostly correspond to the peaks of the initial random gaussian field x- Figure Q] 
shows a snapshot of <j) and x fluctuations shortly after the <j) fluctuations have started to grow. The fluctuations of 
(j) are for the most part in the same places as the fluctuations of %; the correlation between S(j) 2 and 8x 2 in this plot 
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FIG. 3: Values of the field \ on 



a ID slice through the lattice. 
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FIG. 4: Values of field fluctuations 5<j> (red, solid) and 5x (blue, dashed) on a ID slice through the lattice. Because the scale 
of S\ is much larger than Scj) at the time shown, fluctuations of both fields are normalized to the same standard deviation. 

is 0.72. However, the oscillation frequencies are different for the two fields so the fluctuations are in general not in 
phase, as can be seen in the plot. 

Finally, we considered the field statistics for these field fluctuations. It is known that fluctuations produced during 
preheating are non-Gaussian (see e.g. [H). We measure the Gaussianity of the fields through the kurtosis, 3(/ 2 ) 2 / (/ 4 ). 
It is a necessary but not sufficient condition for Gaussianity that this quantity be equal to one. Figure [5] shows the 
evolution of this quantity for fields <fi and x- As noted in [4] the fluctuations become non-Gaussian during preheating 
and begin slowly tending towards Gaussianity thereafter. We see here that Gaussianity is restored after a time of 
order of a thousand inflaton oscillations. 




FIG. 5: Kurtosis for the fields <f> (red, solid) and \ (blue, dashed). The plot on the left shows the development of non- Gaussianity 
during preheating. The plot on the right shows a much larger time scale on which the return to Gaussianity can be seen. 

We can understand this timing with a rough analytical estimate. Perturbatively the time scale for rescattering is 
given by t r ~ l/(na) where n and a are the number density and scattering cross sections for the gas of particles. 
The cross section can be estimated as a ~ g 2 /k* 2 , where k* ~ ((^A) 1 / 4 ^ is the typical particle momentum produced 
during preheating. To estimate the number density we can take n ~ p/k* ~ g 2 < <p 2 >< \ 2 > /k*. Putting all of 
this together 

\f\d>ot r ~ ; 9 o~. (7) 

(. 9 2 /A) 5 / 4 ((^/0o) 2 )((x/^o) 2 ) 

Shortly after preheating the variances ((</> / <j>o) 2 ) and ((x/4>o) 2 ) are both approximately 10 -3 , so the rescattering time 
is of order 

VXdoU - -, ttttt a ^ 1000. (8) 

V) (200) 5 / 4 I0- 6 v ; 



In other words the time that we found is required to restore Gaussianity has the expected order of magnitude. 
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In a lattice simulation all classical dynamics are automatically taken into account, so the Gaussianity of the fields is 
irrelevant to the accuracy of the simulation. While lattice simulations are excellent for describing preheating and the 
evolution shortly afterwards, they can not be extended to time scales long enough to describe the stages of turbulence 
and thermalization. Many approximation techniques are thus employed to describe these epochs (see e.g. Q and 
references therein), and for these it is often important to know the field statistics. Our results suggest that techniques 
that assume Gaussianity should not be employed immediately after preheating, but they also suggest that it can be 
possible to extend lattice simulations long enough to get through the non-Gaussian stage and thus overlap with the 
subsequent period that can be well described in other ways. 

We would like to thank Lev Kofman for useful discussions and advice. This work was supported by NSF grant 
PHY-0456631. 
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